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We propose a new approach to compute correlators of quantum fields in de Sitter space. It is 
based on nonequilibrium field theory techniques, and exploits de Sitter symmetries so as to partially 
reduce the number of independent variables of n-point functions in a manner that preserves the 
usefulness of a momentum representation, e.g., for writing spatial convolution integrals as simple 
products. In this representation, the two-point function of a scalar field only depends on two 
physical momenta, and the corresponding Schwinger-Dyson evolution equations take the form of 
momentum flow equations. Moreover, standard diagrammatic rules can be entirely formulated in 
this representation. The method is suitable for analytical approximations as well as numerical 
implementations. In forthcoming publications, we apply it to resum infrared logarithmic terms 
appearing in the perturbative calculation of vertex and correlation functions. 
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I. INTRODUCTION 



The study of quantum fields in de Sitter space is a 
topic of timely interest. Although the issue of comput- 
ing radiative corrections in curved spaces in general and 
in de Sitter space in particular is a rather old topic (see 
e.g. Q), it has seen a renewed interest in the last decade 
with strong motivations from recent cosmological obser- 
vations. In particular, the impressive success of the in- 
fiationary paradigm in the early Universe 0, Q and the 
observation of the recent acceleration of the Universe [1] 
motivate one to better understand quantum field theory 
(QFT) in expanding space-times. 

A fundamental question is the so-called trans- 
Planckian issue [5], i.e., the issue of an effective decou- 
pling between high- and low-energy physics, which is at 
the root of the concept of effective QFT. If decoupling is 
rather well understood in flat space-time [1, 0] , the sit- 
uation is much less clear in expanding universes, where 
gravitational redshift induces a kinematical correlation 
between infrared (IR) and (arbitrarily high) ultraviolet 
(UV) modes, for which an effective description in terms 
of a fixed background geometry may not be appropri- 
ate. Clear light is shed on this issue when considering a 
lattice formulation of QFT in expanding universes 
[To| . When working with a fixed number of lattice points, 
their density decreases as the universe expands, thereby 
limiting the number of e-foldings of the simulation dur- 



ing which some useful information can be extracted^ ; see, 
e-g-, To avoid this dilution, one is led to consider a 
lattice where the number of sites increases as expansion 
proceeds. It is then a practical question how to initialize 
these incoming degrees of freedom !J] and how to couple 
them with preexistent configurations. 

Besides this issue, one should also consider nonper- 
turbative techniques. Radiative corrections in cosmo- 
logical Friedmann- Robertson- Walker (FRW) space-times 
have been addressed in a variety of field theories d, [Tsl - 

f, mainly based on the perturbative loop expansion 
. Then loop diagrams typically exhibit secular terms, 
which grow as powers of the number of e-folds, and which 
can turn into severe infrared diver gen ces when the field 
is light in units of the Hubble rate (l4 [igI [iTj . Spurious 
secular terms are typical of perturbative approaches for 
nonequilibrium (time dependent) problems in QFT for- 
mulated in flat space-time [l^ . They prevent the study of 
late time evolution and must be resummed to get mean- 
ingful results. Similarly IR divergences naturally arise in 
situations with light bosonic degrees of freedom, such as 
scalar or gauge fields at high temperatures [l^ or near a 
second order phase transition [tI]. They usually signal a 
deficiency of the perturbative approach^. This calls for 
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^ In addition, when working with a comoving lattice, the bare 
parameters of the theory must depend on the cosmological time 
in order to keep the renormalized parameters constant at a given 
physical scale [TUl . 

^ Note, however, that in the context of inflationary cosmology, the 
authors of Refs. |20|| have argued that IR divergences do not 
affect gauge-invariant observables; see also [2l{ . 
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resummation techniques and/or nonperturbative meth- 
ods possibly involving numerical techniques. 

A number of methods have been developed over the 
years to deal with these issues, such as renormalization 
group dill or two-particle-irreducible (2PI) 
techniques. Since the cosmological context is a nonequi- 
librium setup, the most appropriate tools are those of 
nonequilibrium QFT. Techniques such as the dynamical 
renormalization group 24] , or the 2PI formalism [11 [H- 
[It'I can be formulated for expanding space-times. Intro- 
ducing conformal time and comoving spatial coordinates 
as well as conformally rescaled fields, the relevant equa- 
tions for an interacting scalar field theory for instance ac- 
tually very much resemble their Minkowski counterparts. 
The expansion is only manifest in an additional time- 
dependent mass term invol ving both the expansion rate 
and acceleration, see e.g. [ij]. This apparently allows 
one to use the (numerical) tools developed for nonequi- 
librium QFT in fiat space-time directly in this context. 
However, as discussed above, the gravitational redshift 
actually limits the numerical simulations to a low number 
of c- folds [12]. Therefore one must look for approaches 
that avoid this limitation. 

Because of its larger degree of symmetry and of its rele- 
vance to inflationary cosmology, de Sitter space-time has 
been much investigated [Tsl. [2ll [28l - l36j . In this context, 
a very efficient effective description for the nonpertur- 
bative dynamics of IR modes has been devised, the so- 
called stochastic approach 16], which has been shown to 
actually resum the leading IR logarithms of perturbation 
theory to all orders [l3|. It is, however, desirable to go 
beyond this effective description for various reasons, e.g. 
for computing corrections to the stochastic approach, or 
to address specific issues outside its domain of applicabil- 
ity, such as e.g. decoupling, which requires a dynamical 
description of both IR and UV modes. An interesting 
proposal to systematically include perturbative correc- 
tions to the results of the stochastic approach in the con- 
text of Euclidean de Sitter space has been put forward in 
[Ulli]. For scalar field theories, the large- TV |lll,|39|. or 
Hartree [l^, |4l| approximations provide nonperturbative 
approaches capable of taking into account the full cou- 
pled dynamics of IR and UV modes. These essentially 
amount to a local mass resummation and are well suited 
to describe dynamical mass generation in de Sitter. Go- 
ing beyond such mean-field-like descriptions typically re- 
quires one to treat nonlocal integro-differential equations 
[l3 |. Standard nonequilibrium techniques to deal with 
the latter suffer the same drawbacks described above for 
generic FRW geometries. 

In the de Sitter geometry, one may hope to use the 
large symmetry group to further simplify the formula- 
tion and overcome the problems mentioned above. For 
instance, a two-point function in comoving momentum 
space in a FRW geometry with flat spatial sections de- 
pends on three variables: two (conformal) times and the 
modulus of a comoving momentum. One thus effectively 
deals with a (l-l-l)-dimensional problem 12]. In de Sitter 



space, two-point correlators in real space only depend on 
one variable, the de Sitter invariant distance. However, 
the equations of motion for correlators, which typically 
involve nonlocal convolution integrals are difficult to for- 
mulate in real space in a way suitable for both analyti- 
cal and numerical calculations. Instead these are conve- 
niently formulated in Fourier space where, however, the 
full de Sitter symmetry is not transparent and difficult 
to exploit. 

Here, we propose an intermediate approach where we 
exploit only partially the full de Sitter group. Our 
approach is a momentum representation in that con- 
volutions and loop integrals keep a simple form, suit- 
able for numerical implementations as well as for sim- 
plified analytic treatments. It generalizes the so-called 
p-representation introduced in |42| , and exploits the fact 
that the expanding de Sitter space-time can be equiva- 
lently seen as time dependent and spatially homogeneous, 
or as stationary but spatially inhomogeneous. The sub- 
group which combines these two symmetries implies that 
the two-point correlation function (or the vertex func- 
tion) of a scalar field can be expressed in terms of only 
two physical momenta. We then show that the relevant 
Schwinger-Dyson (SD) equations for two-point functions 
take a particularly simple form in terms of these vari- 
ables: the time evolution equation formally becomes flow 
equations in physical momentum space. In addition, this 
effectively reduces to a (0 -I- 1) -dimensional problem. 

We thus see that this approach combines mathematical 
simplifications with physical insight. We believe it pro- 
vides a solution to some aspects of the trans-Planckian 
issue in de Sitter since it allows, in particular, for a nu- 
merical implementation of the basic equations of QFT on 
a grid in physical momentum with no need for adding new 
degrees of freedom as expansion proceeds. Furthermore, 
it opens the possibility of performing numerical calcula- 
tions without being limited by the number of e-folds and 
thus offers a possible way to study various nonperturba- 
tive issues in de Sitter space. Let us finally mention that 
the p-representation is also of relevance in the context of 
(analog) black- hole physics and Hawking radiation [43j . 
where it naturally arises in discussing dispersive and dis- 
sipative effects in Lorentz violating theories [H |3] ■ 

We discuss in detail the physical momentum repre- 
sentation of the original QFT equations in the in-in, 
or closed-time-path formalism in Sec. |TT1 We show, in 
Sec. mil how standard diagrammatic rules for the cal- 
culation of two-point vertex functions can be systemati- 
cally formulated in the p-representation. We then illus- 
trate its usefulness when adopting various approximation 
schemes, such as the perturbative loop expansion and the 
nonperturbative 1/iV-expansion, discussed in Sees. IIVI 
and |V] respectively. Finally, we point out in Sec. IVII 
that the approach is particularly suited for nonpertur- 
bative (resummed) approximation schemes based on the 
2PI formalism. Additional material concerning the re- 
formulation in the p-representation of the in-in closed 
contour integral, of higher order correlation and vertex 
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functions, and of the auxiliary field formulation of the 
expansion are presented in the Appendices. 



II. p-REPRESENTATION OF DE SITTER 
CORRELATORS 



We consider a scalar field ip{x) in the expanding 
Poincare patch of de Sitter space with D — d + I dimen- 
sions and with Hubble scale H — 1. The line element is 
given by 



= -(1 - x^)dt^ - 2x • rfxd< + rfx^ 



(1) 



In the first line, we used the conformal time 77 and co- 
moving spatial coordinates X, referred to as comoving 
coordinates in the following. In the second line, we 
used the cosmological time t and the Lemaitre-Painleve- 
Gullstrand (or physical) spatial coordinates x, hereafter 
named PG coordinates. These coordinate systems are 



with 



related through x — a{rf)~K. and a{rj) 
t e M. 

The comoving coordinates exhibit the homogeneous 
and expanding character of de Sitter space, whereas the 
PG coordinates establish that it is also stationary but 
inhomogeneous. These two facets of de Sitter space are 
at the very origin of the p-representation. References 
(4^ HI] present a detailed discussion of the group theo- 
retical foundation of the latter, which is related to the 
affine subgroup of the de Sitter group'^. 



A. Two-point correlators 

Just like general nonequilibrium quantum systems, 
quantum field theories on FRW geometries are most 
conveniently formulated in the so-called Schwinger- 
Keldysh — also dubbed in-in — formalism [ISj, 2^, 45-48] , 
that is on a closed contour C in the time coordinate; see 
e-g- [13) m, nil- The appropriate contour for conformal 
time is depicted in Fig. [1] and is discussed in Appendix 
El The various components of n-point correlators are de- 
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FIG. 1: The closed path C — C U C in conformal time 

= ri. The forward (upper) branch C goes from —00 to 

0~ and the backward (lower) branch C~ goes back from 0~ 

to — CX3. 



^ We warn the interested reader about the different notations used 
here and in |42l |44|| : here, we use lower (upper) case letters for 
physical (como ving ) variables, which is the opposite convention 
of that used in (UQ • 



scribed by means of time-ordered products of field oper- 
ators on the contour. For instance the two-point func- 
tion G{x,x') = {Tc(p{x)(p{x')), where Tc denotes time 
ordering along the contour C, encodes both the statisti- 
cal and spectral correlators'' F{x,x') = ^{{ip{x), (p{x')}) 
and p{x, x') — i{[ip{x), ip{x')]): 

Gix, x') = Fix, x') - '-signcix° - x'°)p(x, x') , (2) 

where the sign function is to be understood on the con- 
tour C; see Appendix El In the rest of this subsection, we 
consider the statistical function F. Everything we write 
equally applies to the spectral function p 

De Sitter invariance ensures that F{x, x') only depend 
on the invariant distance z{x, x'). In the comoving coor- 
dinate system x = (ry,X), 



z{x, x') 



if + if -jX-X' 



(3) 



When using these coordinates, it proves convenient to in- 
troduce conformally rescaled quantities, such as the field 
4>ix) = {vi)ip{x) and its correlators. One has for in- 
stance 



F{x,x') = [a{r,)a{r,')\ 



F,(r?,77',|X-X'|) (4) 



where Fc(77, 77', |X — X'l) — ^({^(a;), ^(a;')}). Introducing 
spatial comoving momentum variables, one writes, with 
the notation /k = / -^2^, 



iKX+iK'X' 



Fe(r/,77',K,K') 



Fe(7,,77MX-X'|j= / e 
Jk,k 

(5) 

Exploiting spatial homogeneity, one gets 

F,(7;,77',K,K') - (2^)'^<5('')(K + K')/'c(?7,ry',if), (6) 
where 

FM,K)^ f d^5e-*K-«i^,(77,7/,|S|). (7) 



Eq. ^ simply expresses the conservation of comoving 
momentum and is valid in any FRW geometry with flat 
spatial sections. 

Let us now see how this relates to PG coordinates. The 
invariant distance reads 



z{x, x') ^ cosh(At) - i (e^^Atx _ e^^*x'^ ^ 



(8) 



where At — t — t' . We conclude that the two-point cor- 
relator can be written as a function of two variables 

Fix,x') = Fp(Ai, |e-^^*x-e^^*x'|) 

e^P"+^P'"Fp(Ai,p,p'), (9) 

p,p' 



Here, {A, B} = AB + BA and [A, B] = AB ~ BA. 
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where we introduced physical momentum variables in the 
second line. Exploiting the fact that the dependence of 
the two-point function on the spatial coordinates is only 
through the combination |e~2'^*x — e2'^*x'|, one easily 
concludes that the Fourier transform reads 



fp(Ai,p,p') 
with 



At 



FP{P,P') 
(10) 



d^se 



P-^Fp(At,|s|). (11) 



Eq. pUj) expresses the conservation and redshift of phys- 
ical momentum, which is nothing but the conservation of 
comoving momentum K = pe* = — p'e* = — K'. Here, 
we have used the fact that the integral in ()11|) is clearly 
a function of At and p — |p|, which we can trade for p 
and p' — pe^*". 

Combining the two equivalent representations Q and 
([2]) of F, we conclude that 

Fp(Ai,p,p') = [0(77)0(77')]^ F,(r/,?7',a(r,)p,a(r,')p'), 

(12) 

and 

FP{P.P') - HTi)a(if'^ F,{T^,f^\a{f^)p). (13) 
Introducing 



F{p,p') = y^Fp{p,p'), 



(14) 



we get the p-representation of the two-point function [42l, 
El 



F,{r,,^',K) = ^F{p,p') 



(15) 



with p = —Krj and p' — Kr]', the physical momenta 
at times 77 and rj' respectively. This relation expresses 
the fact that the comoving representation of de Sit- 
ter correlators has the scaling property Fc{rj,r]\K) — 
Fc{Kr], Kr]', 1)/K. As mentioned previously Eq. p5|) ap- 
plies to the spectral function p as well. For what concerns 
the calculation of two-point correlators, this effectively 
reduces the number of independent variables from three 
to two. 

Finally, we notice that in going from the comoving rep- 
resentation to the p-representation, the closed contour C 
in conformal time is turned into a closed contour C in 
physical momentum, as depicted on Fig. [2] and discussed 
in Appendix [X] In fact one can grab the statistical and 
spectral component of the two-point function in the fol- 
lowing propagator on the momentum contour: 

^ ^ I 

G{p,p') = F{p,p ) - -sign^ip -p')p{p,p') (16) 

where, as before, the sign function is to be understood 
on the contour; see Appendix [XI We have, in particular, 
signc(?y - v') = sign^~(p - p'). 
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FIG. 2: The closed path C = C+ U in the momentum 
variable p — —Kr). The upper branch goes from -l-oo to 
0^ and the lower branch C~ goes back from 0"*" to -|-oo. 



B. Schwinger-Dyson equations 

We now exploit the above considerations to rewrite 
the SD equations for the two-point functions in the p- 
representation. We first define the covariant inverse prop- 
agator on the closed time contour as 



G-^{x,z)G{z,x')^5^^\x,x') 



(17) 



with = J dP zyj ~g{z) — dz^ J d'^z^J ~g{z), where 
the time integral runs along the contour C and with 



<5(^)(x-y) (5c(xO -?/0)(5('^)(x-y) 



(18) 

the covariant Dirac distribution on the contour, defined 
as J^S^^'>{x,z)f{z) = f{x), see Appendix El 

Schwinger-Dyson equations are obtained by introduc- 
ing the covariant self-energy as 



G ^{x,x) — Gq^{x,x') — Y,{x,x') 



(19) 



where the covariant free inverse propagator is given by 
the quadratic part of the classical action S[ip]: 



zCq (3^, X ) 



where 



S(p{x)Sip{x') 



Sip{x) y^^^6ip{x) 



{a,-mls)S^-^\x,x'), 
(20) 



(21) 



defines a covariant functional derivative. To fix the ideas 
we choose here, in the second equality, the standard form 
of the inverse propagator of a scalar field with standard 
kinetic and mass terms. Here, 



1 



is the covariant Laplace operator and 



Ms 



^R = m'^ + d{d + 1)^ 



(22) 



(23) 



is the effective square mass with m the tree-level mass 
and ^ the coupling to curvature R — d{d + 1). 
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Extracting a possible local contribution to the self- 
energy^ (note that de Sitter symmetry imposes that the 
local term a be constant) , 

Y.{x,x') = -ia5'^'^\x,x') + T.,^i{x,x'), (24) 

the SD equation reads 



[U^- M'^]G{x,x') = i5^^\x,x')+i I T.^i{x,z)G{z,x'), 

(25) 



where we defined 



(26) 



1. Comoving representation 



Let us now specify to the comoving representation. As 
before, we define the inverse comoving propagator with 
conformal rescaling factors: 

G-\x,x') = [a{^)a{v')]-'^ G^\v,v\ |X - X'|) (27) 
such that, in comoving momentum space, Eq. (|17p reads 

Gz\ti, e, K)G,{^, 77', K) ^Sciv- v') (28) 



The conformally rescaled self-energy is defined accord- 
ingly: 

S(x, x') = HrMv'T"^ ^civ, V', |X - X'l). (29) 
and Eq. ()24p becomes, in comoving momentum space, 
Ec(77, V',K)^ -iaa\ij)6c[il - v') + ^fiv, v' , K). (30) 



where we used —g{x) — a^irj). Finally, Eq. (|25|) takes 
the form 



,2 1 



Gc{v,v\K)^-iSc{v-v') 



•J. 



-I l^d^Ef{r,,^,K)Gc{i,ij\K), 



where 



(31) 



(32) 



^ We extract a local term to take into account possible local (tad- 
pole diagrams) contributions to the self-energy. It is to be em- 
phasized that the local contribution to the self-energy is modi- 
fied by the renormalization of UV divergences. In particular, in 
D = A, one expects an additional contribution ~ DxS^*"' {x, x') 
corresponding to field strength renormalization [2^1 . 



In order to write the explicit form of the time integrals 
along the closed contour, we use the standard decompo- 
sition of (nonlocal) two-point functions [l^ 



Gciv, v', K) = /'c(?7, 77', K) - -signc(?7 - 77')Pc('7, K) 



and 



(33) 



i:f{r^, rl, K) - (77, r;', K) - -signc(ry - v'WM v' , K). 

_ (34) 

It is a straightforward exercise to show that [1^ 



d^A{Ti,^)B{£„r^') = 

r diApirj,OBF{C,v')~ r dCAFirj,OBp{^,rj') 

J —00 'J —00 



signc(r;-?7') / d^ Ap{rj,OBp{^,rj'), 



(35) 



so that the SD equations on the time contour read [1/ 
.,2 1 1 



d' + K' ^ 



d^tfiv,^,K)p,{^,T^',K) 
dCmv,^,K)F,i^,r^',K),{36) 



,2 1 



di1:P{7^,^,K)p,{tri'.K). (37) 



These are nonlinear integro-differential equations. They 
are nonlocal and causal since they involve memory inte- 
grals over the whole past history of the system. 

Being second order in time, these equations must be 
supplemented by initial data for the functions Fc and 
pc and their first two derivatives e.g. at 77 = 77' — > 
—00. The initial conditions for the spectral function are 
given by equal-time commutation relations: pc{ri,r], K) = 
dA'Pc{v,v',K)\n=r,' = and ^)|r)=r,' = 1- 

The statistical function contains the information about 
the (quantum) state of the system. Renormalizability (or 
Hadamard conditions) select the so-called Bunch-Davies 
vacuum state as the only viable de Sitter invariant state. 
In the infinite past (subhorizon limit), the latter reduces 
to the corresponding Minkowsky vacuum state of the in- 
teracting theory, which can still be a quite complicated 
state. However, if we define the interacting theory by 
means of an adiabatic switch on of the interaction from 
the infinite remote past, we may identify the state at 
r] — > —00 as the free Bunch-Davies vacuum, character- 
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ized by 



r]—ri — — oo 



'q—iy — — oo 



rj— fj' — ¥ — 00 



1 

0, 
K 



(38) 



As discussed in the introduction, the SD equations 
in the comoving representation, Eqs. (|36|) (|37|) . have a 
very similar structure as nonequiUbrium evolution equa- 
tions for a scalar field in flat space-time: the only place 
where expansion enters is the time-dependent mass term 
oc 1/rf'. In this representation, the calculation of de Sit- 
ter correlators is fully expressed as an initial value prob- 
lem. Exploiting the spatial homogenetity and isotropy 
of de Sitter geometry in comoving coordinates, the prob- 
lem effectively reduces to the calculation of a (1 + 1)- 
dimensional nonequilibrium two-point correlator. 

We can readily see the difficulties mentioned in the 
introduction on Eqs. (I55|) - (I571) . If one chooses a dis- 
cretization on a (fixed) grid in comoving coordinates, the 
growing mass term eventually becomes larger than the 
comoving momentum cutoff. It thus requires an initially 
extremely fine lattice in order to resolve the inverse mass 
during a large number of e- folds. At the same time one 
wants an as large as possible spatial volume in order to 
correctly describe IR physics. In practice, this approach 
is bounded to only a few e- folds fl3 |. 

Moreover, as already emphasized, such a discretization 
does not have a continuum limit in Z? = 4. A more appro- 
priate choice in this respect is to discretize the system in 
proper physical coordinates. However, in that case, the 
number of comoving modes K involved in the simulation 
increases with time. The will to correctly describe IR 
physics (which requires a large volume) for a long time 
(which eventually requires a large number of degrees of 
freedom) leads to a similar difficulty as in the previous 
case. Another issue in that case is that one needs to 
supplement the evolution equations with an ad hoc spec- 
ification of how to initialize the new degrees of freedom 
which constantly enter the system. 

Let us now discuss how SD equations can be formu- 
lated in the p-representation and show how this solves 
the above issues. 



2. p -representation 

As already mentioned the contour C in conformal time 
can be traded for a contour C in momentum, as depicted 
in Fig. [5] above. We first define the p-representation of 
the inverse propagator as 



with p = —Ki] and p' = —Krj', such that Eq. (l28l) be- 
comes 

j^dsG-\p,s)G{s,p')=5^{p-p'), (40) 

where 5c{fl) is the delta function on C, defined such that 
j^ds5^{p — s)J{s) = f{p), see Appendix Rl Notice, in 
particular, the relation Sc{r] — ij') — —KS^{p — p')- 

Assuming that the self-energy scales as the inverse 
propagator^. 



t,{r,,r^\K)^K^t{p,p'), 



(41) 



one can write 



Y.{p,p ) = ta-^ 



+ tnl{p,p'), (42) 



where the nonlocal contribution is defined as 
tf{T],r]',K) = K^t^i{p,p'), and the SD equation 
can be rewritten fully in the p-representation as 



5; + 1 - 



G{p,p') = i6f^{p-p') 

+ iJ^ds±ni{p,s)G{s,p'),i4S) 



As before, the contour integral can be written explicitly. 
We write 

G{p,p') = F{p,p') - ^sign^ip -p')p{p,p'), (44) 

^ ^ 1 

Sni(p,p') = T.f{p,p') - -signg"(p -p')I]p(p,p'), (45) 

and similarly for any nonlocal two-point function on the 
contour C. Eq. (|35|) becomes 



i / dsA{p, s)B{s,p') = 



poo poo 

/ ds Ap{p,s)Bf{s,p') - dsAF{p,s)Bp{s,p') 



sign^ip^p') I dsAp{p,s)Bp{s,p'), (46) 



G-\v,fj\K)=K'G-'(p,p'), 



(39) 



® There is a freedom in the choice of the p-representation of 
the inverse propagator. For instance, in Ref. [44| . the co- 
variant inverse propagator and propagator are treated on an 
equal footing. In the present notations, the authors of [43 ] 
thus introduce the following function, sec Eq. (|4]l: E(x, x') = 

_ d-l 

[a{'q)a{'q')] 2 Sabp(^i^')|X — X'|), which admits the p- 
representation, see Eq. Ijl5|l Sabp(^i J?'; ^) = ^ABp(Pi p')/^- 
With this choice, the convolutions in the physical momentum 
variable, e.g. in Eqs. I I47I1 - II48I I. involve a nontrivial measure 
J ds/s^. Here, we treat the covariant inverse propagators and 
propagators differently. Our choice is such that convolutions in 
the physical momenta involve a trivial measure. Our self-energy 
is related to that of Q by Sabp(p,p') = (pp')^^(P.P')- 
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and the SD equations thus read 

,,2 1 1 



d 



p 



1 - 



F{p,p') = / dsT.F{p,s)p{s,p') 



ds Sp(p,s)F(s,p')> 



(47) 

ds tp{p,s)p{s,p'). 

(48) 



The "initial" data are to be specified at p — p' —> 
oo. Commutation relations imply p(p,p')lp=p — 
dpdp'p{p,p')\p=p' = and dpp{p,p')\p=p' = -1 for the 
spectral function and the choice of the free Bunch-Davies 
vacuum at large momentum — keeping in mind an adia- 
batic switching on of the interaction — means 



F{p,p') 
dpFip,p') 
dpdp'F{p,p') 



p—p'—^OQ 
p—p'^OG 



1 

2 



1 

2 



(49) 



for the statistical correlator. 

Eqs. (|47)) (|48l) generalize the evolution equations of 
a free field in the p-representation introduced in E^. 
They provide an important simplification as compared to 
the comoving formulation in that the problem is reduced 
to an effective (0 + l)-dimensional, quantum-mechanical- 
like problem, with two-point correlators depending only 
on two momentum/time variables, instead of an effective 
(1 -|- l)-dimensional problem with two time an one mo- 
mentum variables in the comoving representation. This 
is of particular importance for numerical investigations 
of SD equations on de Sitter space using nonequilibrium 
techniques. 

It is remarkable that in the p-representation the time 
evolution is replaced by a momentum evolution and, as 
a consequence, for given self-energy kernels E_F^p, the SD 
equations turn into integro-differential flow equations in 
physical momentum with a second order derivative. This 
is for instance the case in ordinary perturbation theory, 
where the self-energies at a given order are given func- 
tions of the free propagators. In such a case, the cal- 
culation of the propagator for a given momentum only 
involves higher momenta and Eqs. (|T7|) describe 
how the integration of higher momenta builds up lower 
momenta correlators. However, this is no longer the case 
for self-consistent approximation schemes such as those 
based on the 2PI formalism, where self-energies involve 
momentum (loop) integrals of the full propagators them- 
selves and thus depend on the latter at all momenta. In 
such a case one can envisage obtaining a nonperturbative 
solution by means of iterative techniques [49| . 



III. DIAGRAMMATIC RULES 

The considerations of the previous section rely on 
the scaling assumption (j4ip . which ensures that the p- 
representation closes. However the self-energy is gener- 
ated by the field (self-)interactions and one has to check 
that this scaling relation is actually satisfied whenever 
the scaling ([T5|) for the correlator is true. To do so, in the 
next subsection, we shall analyze the diagrammatic rep- 
resentation of the self-energy and derive Feynman rules 
in the p-representation. The case of higher order correla- 
tion or vertex functions is briefiy discussed in Appendix 
IB] Finally we give some explicit examples of self-energies 
in the p-representation in both perturbative and nonper- 
turbative approximation schemes. 



A. Two-point functions 

We start from the usual diagrammatic rules for the 
self-energy E(a;,a;') in the covariant formulation. We are 
thus concerned with connected one-particle-irreducible 
(IPI) diagrams with two external (amputated i.e. not 
associated to a propagator) legs. We do not need to 
specify any particular interaction term. We assume 
nonderivative — but otherwise arbitrary — polynomial in- 
teractions. Consider a given IPI diagram: each inter- 
nal line contributes a G{zi,z'j) with zi and z'^ the end- 
points of the ith. line; each vertex with n legs contributes 

^ rife^i <5^'°H-^fc) ^k+i) "^'^^^ ^(j ■ ■ ■ T^ii the space-time co- 
ordinates associated with each leg of the jth vertex; fi- 
nally all coordinates associated with vertices must be in- 
tegrated over with the covariant measure, J^, except for 
the two coordinates x and x' associated with the two ex- 
ternal legs. For a diagram with / internal lines there are 
21 space-time coordinates to be integrated over. 

Let us now move on to the comoving representation 
in momentum space, i.e. to the diagrammatic rules 
for Yjc{rj,ri' ,K). First, there is now an overall factor 
[0(77)0(77')] T~ = (7777') i~ from the definition of the 
conformally rescaled self-energy, Eq. (|29|) . To each inter- 
nal line of the diagram under consideration is associated 
a comoving momentum Qi, to be integrated over, and 
two conformal time endpoints "qi and 77^. Each such line 

contributes a factor {r]iri'^)^~Gc{'r]i:Vi:Qi)- Each vertex 
with n legs contributes a comoving momentum conserva- 
tion factor (27r)'*(5('*) iJ27=i Qi), where the sum runs over 
all momenta entering the vertex. One of these Dirac 
factors ensures the total comoving momentum conserva- 
tion, as e.g. in Eq. (|6]), and is extracted in the definition 
of Ec. Therefore, for a diagram with V vertices, there 
are V — I momentum conservation factors. The number 
of independent momentum integrations (loops) is thus 
L — I — V + 1, the usual relation. Each vertex not at- 
tached to an external endpoint contributes an integral 
over a conformal time variable drjia^ (rji). There are 
V — 2 such vertices for nonlocal contributions to the self- 
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energy (local terms are to be treated separately). 

We next translate the above rules to the p- 
representation. To do so, we rescale all — internal and 
external — conformal time variables as well as internal 
momenta with the external comoving momentum K: 
p — —Krj, p' = —Krj' for time variables associated with 
external legs, pi = —Krji, p'^ = —Kri[ for times variables 
associated with internal vertices, and — K(\i for in- 
ternal momenta. Each line factor then reads 

(50) 

We see that although some cndpoints rji , rjl of some lines 
are actually equal to the external ones 77,77', the line fac- 
tors can be entirely written in terms of p, p' thanks to the 
fact that only ratios, e.g. 77/77.; — p/pi, occur. Now, each 
of the V — 1 comoving momentum conservation terms 



(51) 



contributes a factor K^"^ and each integral over confor- 
mal times associated to internal vertices (not attached to 
an external leg) 



dr]ia°{rii) 



drji 



i-mV 



S (52) 



c P. 



contributes a factor^ K'''. 

A nonlocal contribution to E is such that the two ex- 
ternal legs are not attached to the same vertex. A dia- 
gram with V vertices has thus V — 2 internal vertices and 
contributes a term 



V-2) 



Kd(V-l) 



X (7777' 



X {pp') 



(53) 



times a function of p and p' only. We emphasize that this 
is independent of the number of vertices V and thus of 
the particular diagram under consideration. This demon- 
strates the scaling of self-energy, Eq. (|4T|) . at any or- 
der of perturbation theory as a consequence of the scaling 

m- 

For a local contribution to the self-energy, both exter- 
nal lines are attached to the same vertex and there are 
thus V — 1 internal vertices and an extra 5c{ri — T]')/a^{rj). 
Altogether one gets an overall factor 



as required for the p-representation, see pop and ([^^ 



(54) 



Here, the minus sign on the right-hand side is to be included in 
the diagrammatic rule for integrating over internal vertices. It 
simply reflects the orientation of the momentum contour. 



Finally, we emphasize that the previous analysis shows 
that the diagrammatic rules in the p-representation are 
the same as those in the comoving representation with 
the generic replacements Q — q for all momenta (in- 
cluding the external one for which one has K — e), 
—77 — > p for all time variables and Gc — > G/g for 
all propagator lines, see Eq. ((5(1)) . In the next subsec- 
tions, we give explicit examples for various approxima- 
tion schemes. The diagrammatic rules for higher corre- 
lation and vertex functions in the p-representation are 
discussed in Appendix |B] 



IV. LOOP EXPANSION 

We illustrate the above considerations for an 0{N) the- 
ory with quartic coupling. With the definitions ([^ and 
(1231), the classical action reads 



(55) 



where a = 1, . . . , A'^ and a summation over repeated in- 
dices is understood. We consider the symmetric phase, 
{(fa) = 0, for which the propagator and self-energy are 
diagonal: Gab = SabG and Y^ab = (^afcS. 

As a first example of an approximation scheme, we 
consider the standard loop expansion which, in the case 
under consideration, is equivalent to a coupling expan- 
sion. We write the formal series 

S = + ... (56) 

where E*^"-' ^ ©(A") is the 71-loop order contribution. 



A. One loop 

Let us first recall the result in the covariant formula- 
tion. At one-loop order there is only a local contribution 
given by the tadpole diagram of Fig. [3] 



with 



gGQ{x,x). 



(57) 



(58) 



where we defined g = X{N + 2)/6N and where Go de- 
notes the free (covariant) propagator; see Eq. (^U)) . The 
de Sitter symmetry group guarantees that Go{x,x) only 
depends on the invariant distance z{x,x) = 1. The mass 
shift (7(1) is thus a constant. Applying the diagrammatic 
rules directly in the p-representation, we obtain 



9 / P 
'q 



d^iGoiqp^qp) 



9 



Fo{q,q) 



(59) 



where we made the change of variable q — J> q/p and used 
the representation (jl6p for the free propagator on the 
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contour in the second equality. This is readily seen to 
coincide with the above expression^ Eq. (|58p . using Eqs. 
Cni) and (HI). 




FIG. 3: The one-loop tadpole contribution to the self-energy 
'E{x,x'). The black dot denotes an interaction vertex and 
the line in the loop represents the free propagator Go (a;, a;). 
A similar diagram describes local contributions in the 
expansion. In that case, the line represents the leading order 
(large- A'') propagator G{x,x). 



B. Two loop 



the contour C explicitly using Eq. (l46l) . one also checks 
that crQ^-* is real as expected^: 



The nonlocal contribution reads, in the covariant rep- 
resentation. 



(64) 



where g' ~ —X'^{N + 2)/18A^^, or, in the comoving rep- 
resentation. 



(65) 



JQ,L 

where R — \Ke + Q + L| with e an arbitrary unit vector. 



At two loop there is both a local and a nonlocal contri- 
bution. They are depicted in Figs. |3] and [5] respectively: 

x') = -z(t(2) j(^) {x, x') + [x, x'). (60) 

The local contribution is a first example with an internal 
vertex. It reads, in the p-representation, 

f^(2)=ig2 f ±_ )d-iG|(fcPvfcf) fd-iGoiqsM 

Jc A k'^ Jq q ' 

(61) 

Applying the changes of variables q — ^ q/s, k — i> k/p 




FIG. 4: The local two-loop contribution to E(a;,x'). 
and s ps and using (j59p . it can be rewritten as 

Jc s Jk 

which clearly shows that it is indeed a constant, as re- 
quired by de Sitter symmetry. Writing the integral on 




FIG. 5: The nonlocal two-loop contribution to 'E{x,x'). 

A direct application of the diagrammatic rules in the 
p-representation gives 

f,(2). n nd-sf Go{qp,qp')G(){lp,¥)Go{rp,rp') 



qi 



qlr 



(66) 

where r — |e + q + 1|. Eq. (|66| is easily checked to 
coincide with the previous expressions (|M1) or when 
converted in the appropriate representation. The explicit 

(2) (2) 

expressions of the component T,p and Sp are also easily 
obtained: the product GqGqGo under the integral gives 
rise to the combinations (keeping the same momentum 

arguments) FqFoFq - fFoPoPo for and SFq^oPo - 



V. 1/A EXPANSION 

Let us now consider an example of a nonperturbative 
approximation scheme, the l/N expansion. The latter 
is a powerful tool to describe nontrivial IR physics in 



^ An alternative expression is 

Note that the divergent integral in Eq. H59| l needs to be regulated. 

Note also that a cutoff on comoving momenta would lead to a ^(2) _ „2(;(7'^' 

f f°° ds Foik,s)po(k,s) 

time dependent result and would thus be inconsistent with the Jk Ji^ k 

de Sitter symmetry. 
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situations where perturbation theory fails. Exact results 
for IR de Sitter correlators have been recently obtained 
in the large- iV limit [U [13], which reveals interesting 
phenomena such as radiative symmetry restoration, or 
the generation of anomalous dimensions. We write the 
formal series in 1/A^ 



(67) 



where E^^" - 0{N°), E^lo _ 0{l/N), etc. The di- 
agrammatics of the 1/7V expansion in flat space-time 
is well known [5ll - [55l |. The generalization to arbitrary 
background geometry is straightforward in the covariant 
formulation. We do not recall the derivations here but 
merely state the results and show how they can be writ- 
ten in the p-representation. We follow the notations of 
Ref. [54]- 



A. Leading order 

The leading order (LO) contribution is a simple tad- 
pole diagram; see Fig. [H It is local and has a similar 
structure to the one-loop result discussed above. The es- 
sential difference is that the tadpole loop is given self- 
consistently — hence the nonperturbative nature of the 
approximation scheme — in terms of the full LO propa- 
gator G. One has 



with 



A 



a^^^-Gix,x) 
6 



(68) 



(69) 



The LO propagator G is defined by [see Eq. (1^5)) ] 

{a^ Mlo)G{x,x') ^i6^^\x,x') (70) 

with 

Mlo^rnls + cT^°. (71) 
The p-representation of the LO approximation reads 



LO 



Fiq,q) 



(72) 



The diagrammatic 1/A^ expansion can be expressed 
fully in terms of the LO propagator G, which resums 
the infinite series of so-called daisy and superdaisy tad- 
pole diagrams [s^ . Alternatively, it proves convenient to 
introduce an auxiliary composite field y cx ipa^a in order 
to organize the 1/A^ expansion |5ll-l55|. We shall follow 
the first approach here. The auxiliary field formulation 
is briefly discussed in Appendix [Cl 



B. Next-to- leading order 

The next-to-leading order (NLO) contribution contains 
both a local and a nonlocal part: 

ENL0(^^^') ^ -jaNLO^P)(^^^') + sNLO(^^^'). (73) 
The local part is simply given by 



NLO _ J^rf-r t) - A^LO 



(74) 



and the nonlocal part resums the infinite series of dia- 
grams shown in Figs. [6]and[7l It can be written as 

i:^t''{x,x')^^G{x,x')I{x,x'), (75) 

where the function / resums the infinite series of bubble 
diagrams shown in Fig. [5] through the following integral 
equation 

I{x, x') = n{x, x')+i / n(a;, z)/(z, x'), (76) 



(77) 



x' + ... 



with the elementary one-loop bubble 



A 



Uix,x')^-'jG^x,x') 




.r' + X 




x! + X 




FIG. 6: The infinite series of bubble diagrams contributing 
to the function I{x,x'), Eq. (|76p . The black dots correspond 
to interaction vertices whereas the crosses denote the end- 
points of the function. The elementary bubble is given by the 
function 11(2;, x'), Eq. (|77[1 . Each additional bubble involves a 
summation of field components and thus comes with a factor 
A^, which is compensated by a 1/A'' from the corresponding 
additional vertex. All such diagrams are thus of the same 
order in 1/N. 




FIG. 7; A typical multiloop diagram contributing to the self- 
energy 'E{x,x') at NLO in the 1/N expansion. The latter 
actually resums all diagrams of similar topology with an ar- 
bitrary number of bubbles in the upper part, as described by 
the function I{x,x'); see Fig. [B] 

The NLO contribution can be expressed in the co- 
moving representation by introducing the conformally 
rescaled quantities 

U{x, x') = [a{rj)a{rj')]-'^n,{rj, rj' , |X - X'|) (78) 
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and similarly for /. One obtains, in comoving momentum 
space, 

Ylciv, v', K) = -^ (W)"^ / Ge iv, v\ Q) Gc (^, v', R) , 

(79) 

where R = | Ke + Q | , and 



Iciv,v',K)=Uc{7],r/,K)+i / diUcirj,^,K)Ici^,rj',K). 

Jc 

(80) 

Finally the nonlocal part of the NLO self-energy reads 
SNLO(^, ^'^ K)^^ (W)^ / Ge (77, v', Q) ic {v, v', R)- 



3N 



Q 



(81) 

Using the methods described in previous sections, it is 
easy to check the scaling relations 

Ilc{v,v',K) = Ktl{p,p') , ic{v,v',K) = Ki{p,p') 

(82) 

where p = —Krj and p' = —Kr]', which can be used to 
convert the above equations to the p-representation. One 
gets 

niP,p') - -I iPP')'-^ I (83) 

6 Jq 9 r 

where r = |e + q|. and the function / satisfies the integral 
equation 

i{p,p') ^Il{p,p') -t J^dsll{p,s)i{s,p'). (84) 

Finally the nonlocal part of the NLO self-energy reads 

t^nP,P') = ^ iPP')"^ j \ g {qp, qp') i {rp, rp'). 

(85) 

Again one can check that Eqs. ([55| ((55|) can be obtained 
by direct application of the diagrammatic rules in the 
p-representation as described in the previous section. 

Let us finally write the above equation in terms of the 
explicit components on the momentum contour C. The 
product & in the expression of the function 11 gives 
rise to the combinations FF — jpp for lip and 2Fp for 
Hp. Similarly, the product GI in (1851) gives FIp — jpip 
for E^L° and Pip + pip for Finally, the contour 

integrals in (IMl) are obtained from Eq. (^5)) as 

/>oo 

iF{p,p') - IIf{p,p')~ dsUF{p,s)ip{s,p') 



(86) 



ds flp{p, s)/f(s,p'), 



ip{p,p') = tlp{p,p)+ / dstlp{p,s)Ip{.s,p'). {87) 



We end this section by mentioning that the infinite se- 
ries of bubble diagrams discussed here is actually related 



to the four-point vertex function in the large-iV limit. 
The latter is studied in Ref . [13] , where the above integral 
equations are solved exactly in the limit of IR momenta 
p,p' <C 1, making extensive use of the p-representation. 



VI. 2PI APPROXIMATION SCHEMES 

An important class of approximation schemes is based 
on 2PI functional methods [56l - l59l |. These provide sys- 
tematic infinite resummations of selective sets of pertur- 
bative contributions and have proven a very useful tool 
in recent years to resum infrared divergences of bosonic 
theories in fiat space-time at very high temperatures Il9| 
or secular divergences of nonequilibrium field theory |l8| . 
It has been shown that these methods are also useful in 
dealing with infrared and secular issues in de Sitter ge- 
ometry [ill, [2^ [23, H^- We briefly recall the main in- 
gredient of the 2PI formalism in a nonequilibrium setup 
[isl [23I and show how it can be formulated in the p- 
representation. 

2PI self-consistent approximation schemes are based 
on truncations or systematic expansions of the 2PI ef- 
fective action, r[(/), G], a functional of both the one- and 
the two-point correlation functions of the theory in the 
quantum state under consideration 4>a(x) = {ipa{x)) and 
Gab(x,x') = {Tcipa{x)ipb(x')) . It Can be parametrized as 



r[0,G] = + -TrLnG-1 + -TrGo 'G + ri„t[0, G], 



where both the trace Tr and the logarithm Ln are to be 
understood in the functional sense. Here S is the classi- 
cal action, iG^^ is the inverse free covariant propagator 
and Pint can be represented as the infinite sum of closed 
2PI diagrams with lines G and vertices — including two- 
leg vertices — given by the shifted action S[(l) -\- (p\. Such 
diagrams with lines given by the exact propagator of the 
theory instead of the perturbative one are called skeleton 
diagrams. 

All vertex and correlation functions of the theory can 
be obtained from functional derivatives of the 2PI effec- 
tive action evaluated at the solution of the equations of 
motion for both 6 and G: 



5MG] 



0. 



<5,r[0,G] 



0, 



5(t)a{x) SGab{x,x') 

where we define the covariant functional derivatives [2 



(89) 



1 



Sc/yaix) ^J-g{x) Ha{x) ' 

Sr. 1 1 



(90) 



SGabix,x') \/-g{x) ^J-g{x') SGabix,x') 



TV (91) 



In particular, using the parametrization (I88p . the second 
equation gives the SD equation 



G-l = Gn 



(92) 
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with 



with 



Y,abix,x') = 2 



. Squint [0, G] 
5Gba{x',x) 



(93) 



The point here is that the self-energy thus obtained is 
typically a nonlinear functional of the full propagator and 
one has to solve Eq. ((92|) self-consistently. This is where 
the nonperturbative nature of this approximation scheme 
enters. 

We now observe that 2PI approximation schemes are 
formulated in terms of two-point functions and skeleton 
diagrams. Thus, since the full propagator has the correct 
scaling ([T5l) . all the considerations of previous sections 
concerning the SD equations and the diagrammatic rules 
in the p-representation hold. The only modification is 
that the free propagator is replaced by the full one, to be 
determined self-consistently by solving the SD equations. 

For simplicity, let us consider 0{N) symmetric states, 
for which (pa = 0, Gab = SabG and Sab = Sab^- In 
that case, the self-energy E can be obtained from the 
functional Tint evaluated in the symmetric configuration 

(j)a = 0, Gab = SabG aS 



(94) 



N 5Gix',x) 
For instance a 2PI loop expansion at two-loop order gives 

ri„t[0 = 0, G] = / G^x, x)-i^-^ I GHx, y), 

(95) 

with g = \{N + 2)/67V and g' = -\^{N + 2)/im'^ de- 
fined previously. One obtains for the self-consistent self- 
energy 

T.{x, x') = -igG{x, x)6^^'> {x, x') + g'G^{x, x') (96) 

These expressions have the same structure as the stan- 
dard IPI one-loop expressions described in the previous 
section^" and can thus be easily formulated in the p- 
representation. It is sufficient to replace free propagators 
by full ones in the diagrammatic rules. 

Similarly, the 2PI 1/A^ expansion |54J at NLO gives, 
up to an unphysical constant, 

ri„t[<A = 0,G] = / G^{x,x) + ^TrLnl?-!, (97) 



where 



iD^^{x,x') 



3iv r 



S'^^\x,x')-iU{x,x') 



(98) 



In the 2PI loop-expansion, the diagram of Fig. [4] is absent be- 
cause of the 2PI character of the diagrammatic expansion, which 
in fact avoids possible double-counting. The missing two-loop 
contribution is now included in the self-consistent propagator. 



U{x,x') = -^G'ix,x'). 
One easily checks that 

^D{x,x') = -^[s^''Hx,x') +il{x,x') 



(99) 



(100) 



(101) 



The corresponding self-consistent self-energy is given by 
S(x, x') ^ -iaa5^^^ {x, x') + -^G{x, x)I{x, x'), (102) 



where 



I{Xtx') = Il{x^x) + i / n(x, 2)7(2;, x') 



with 



<^o = ^(l + ^)G(a:,a:) 



(103) 



We see again that the structure of the equations is the 
same as those discussed in the previous section, the only 
change being that the LO propagator is now replaced 
by the full one. The p-representation of the 2PI 1/A^ 
expansion is thus obviously obtained. 



VII. CONCLUSIONS 

We have developed a systematic method for comput- 
ing correlation and vertex functions of a scalar field in de 
Sitter space. It exploits both the simplifications due to 
de Sitter symmetries and the power of a momentum rep- 
resentation, e.g., for writing spatial convolution integrals 
as simple products. The method relies on the particular 
way momentum redshift is encoded in de Sitter correla- 
tors, see Eq. (|10p. which implies a one-to-one correspon- 
dence between time and physical momentum and allows 
one to trade one for the other. 

This method is particularly well adapted to describ- 
ing two-point functions — for which it reduces the num- 
ber of independent variables from 3 to 2 as compared 
to the comoving momentum representation — and thus to 
all approximation schemes based on the use of the latter. 
This includes standard expansion schemes in QFT such 
as the loop, or the 1/N expansions, but also resummed 
approximation schemes based on 2PI techniques. For 
what concerns two-point correlators, our approach effec- 
tively reduces the problem to a (0 + l)-dimensional one, 
where physical momentum plays the role of the "time" 
variable. 

We emphasize that the resulting equations are well 
suited to analytical approximations as well as to numer- 
ical implementation. This is particularly important for 
studies of infrared/secular issues in de Sitter, which re- 
quire infinite resummations and thus possibly numerical 
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work, or for studies of trans-Planckian issues which may 
require nonperturbative calculations of unequal time (un- 
equal momentum) correlators, in the same spirit as the 
calculation of damping and thcrmalization effects from 
2PI techniques in flat space-time [18,] . 

Work in these directions has been pursued and the re- 
sults will soon be presented fio!, . In \5d\ , we studied 
the four-point vertex function of an 0{N) scalar field. In 
the large- -/V limit, this vertex is given by an infinite series 
of bubble diagrams, each of which exhibits large IR loga- 
rithms, typical of de Sitter space. Exploiting the method 
presented here, we found by analytical means that the re- 
sulting resummation of IR logarithms leads to a modified 
power law in the deep IR, analogous to the generation of 
an anomalous dimension in critical phenomena. 

Finally, we mention that the present method, since it 
does not exploit all de Sitter symmetries, allows one to 
treat deformations of de Sitter which are compatible with 
the p-representation. In fact the usefulness of this repre- 
sentation was first understood in the context of theories 
where the Lorentz group is violated by dispersive or dis- 
sipative effects occurring in the UV sector [Z^, 113 • In- 
teresting extensions of the present work also include the 
discussion of the p-representation for fields of higher spin 
and/or theories with derivative couplings. 
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Appendix A: Time and momentum contours 



The notion of a closed time contour [45|-|48j provides 
a convenient way to deal with the various components 
of 71-point functions in a general noncquilibrium setup. 
Alternative formulations employ tensor structures, such 
as the so-called ± or Keldysh basis |48[ . 

The closed contour C = UC~ in conformal time can 
be described by a mapping 77(5) of the real axis R onto 
the contour C, such that the negative and positive real 
axis are mapped onto the upper and lower branches 
respectively: 



77(3) G for s e 



(Al) 



with 77(5) a monotonously increasing (decreasing) func- 
tion on the negative (positive) real axis. The ordering 
along C is defined as an ordering along R: 

signc('7 - v') = signc(77(s) - r]{s')) = sign(s - s'). (A2) 

The integral of a function /{rj) = f{ri{s)) along the con- 
tour is defined as 



We define the Dirac delta function on C as 

6c{v^r^')^Sc{v{s)~r^is'))^(^^^ 5(3 - s') (A4) 
in such a way that 

(A5) 

Notice the properties 

^signc(77-?/) = 2<5c(?/-7?') (A6) 



and 



(A7) 



The contour C in momentum is defined in an analogous 
way. We introduce the function {K is a constant) 



pis) ^ -Kri{s), 



(A8) 



which maps the negative and positive real axis R^ onto 
the upper and lower branches respectively: 

p{s)eC^ for seR=F. (A9) 

Here, p{s) is a monotonously decreasing (increasing) 
function on the negative (positive) real axis. The or- 
dering along C, defined as 

signg-(p-p') = sign^{p{s)-p{s')) =sign(s-s'), (AlO) 

is such that 

signj(p - p') = signc(?7 - v')- (AH) 
For a given function 

/(r?(s)) = f{p{s)) (A12) 
we define the integral 

dpf{p) = f ds^fipis)) ^^K [ dr^firj). (A13) 

C JR "'S Jc 

We define a Dirac delta function on C as before: 

<5,-(p-p')^<5,-(p(.)-p(.'))- (J) 'sis-s'), (A14) 
such that 

dpS^{p~p')f{p) = / dsS{s - s) /(p(s')) = fip)- 

JR 

(A15) 

One has the properties 

^sign^(p - p') = 25^{p - p') (A16) 



and 



5r{p~p')^TS{p^p') if P,p'eC±. (A17) 



(A3) 



Finally, notice the relation 

5dTl-ri')^-K5^{p-p'). (A18) 
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Appendix B: Higher order correlators and vertex 
functions 



In this appendix, we discuss the consequences of the 
p-representation ([T5]) of two-point functions for higher 
order correlation and vertex functions, as obtained from 
perturbative diagrams. The connected time-ordered n- 
point correlator is 



G'^"^(a;i, . . . ,Xn) = {Tc(p{xi) ■ ■ ■ ip{xn))c 



(Bl) 



and we define the covariant n-point proper vertex func- 
tions r^") from the field expansion of the IPI effective 
action 

(B2) 

Equivalently 



Sip{xi) ■ ■ ■ 5if{Xn) 



(B3) 



'P=4> 



with the covariant functional derivative ((2T|) . The corre- 
sponding conformally rescaled quantities are given by 

G(")(xi, . . . , a;„) = [a{m) . . . a(?7„)]-^G(")(a;i, . . . , a;„) 

(B4) 

and 

r(")(a;i, . . . ,x„) = [a(ryi) . . . a{nn)]'^V^'\xu ■ • ■ , x^). 

(B5) 

The comoving momentum representation is defined as 
G(")(x,)=/ e'^"-iK^-^^G(")(7;„K,), (B6) 

JKi,...,K„ 

in) 

and similarly for Tc , where we denote collectively Xi = 
xi, . . . , a;„ and similarly for ry's and K's. Spatial homo- 
geneity in comoving coordinates implies that 



only V — n time integrals over internal vertices. We thus 
get the overall factor 



Kdiv-1) ii-m) ■ ■ ■ {-Vn)Y 



K 



d+- 



^ x(pi---p„)" 



(B8) 



which generalizes Eq. (j53p for n ^ 2. Here, we introduced 
the variables pi = —Krji = pi/|Ki|, with pi = — |Ki|r/i 
the physical external momenta. The factor (IBSP is in- 
dependent of V and thus of the diagram one considers. 
We conclude that the n-point vertex function admits the 
scaling (or p-representation) 

f(")(7y„X,K,) = i^'^+^f(")(^,K,). (B9) 

A similar analysis for the n-point correlator ci"^ leads 
to the p-representation 

G(")(77„X,K,) = X'*-^G(")(p„KO. (BIO) 

Eqs. (|B9l) and (jBlOp generalize Eqs. dH]) and (HI]) for 
two-point functions (for which a natural choice for the 
scale K isK = |Ki| = IKal). 



Appendix C: The auxiliary fleld formalism 

The 1/A^ expansion discussed in Sec. |V] can be con- 
veniently formulated by introducing an auxiliary field 
X corresponding to the composite operator ip"^ |5ll455| . 
Consider the action 



2 ^ 3iV 2 

ifa (□ - m^sj + —X 



(CI) 

The equation of motion for the auxiliary field x is a con- 
straint equation: x = ^fa^a/QN. Integrating out the 
field X, one recovers the original action ([55)) . The field x 
is a convenient bookkeeping device to count factors of N 
and thus to organize the 1 /N expansion. 



G(")(r;„K,) - {2^f5^''^ K^K, G(")(,7„if,K,). 



(B7) 

(n) 

and similarly for Tc ■ Spatial isotropy further reduces 
the number of independent momentum variables. Here, 
we extracted an explicit arbitrary comoving momen- 
tum scale K from the set of variables and define 
Ki ~ ICi/K. For instance, one can choose to single out 
one momentum: K = |Ki|, or take the symmetric com- 
bination K = ^/Kl + ...+ Kl, etc. 

(n) 

Now, we consider a given diagram contributing to Tc 
with V internal vertices. We apply the diagrammatic 
rules recalled in Sec. IIIII and perform the same rescaling 
of all — internal and external — variables as in Sec. IIIII 
using the scale K defined above. Since there are now n 
vertices attached to the endpoints Xi, . . . ,Xm there are 



FIG. 8: The NLO contribution to the self-energy T,{x,x') in 
the 1/A'' expansion in the auxiliary field formulation. The 
plain line represents the propagator G of the field ip whereas 
the dashed line denotes the propagator D of the auxiliary 
composite field x ^ 'P^ ■ The latter resums the infinite series 
of bubble diagrams depicted in Fig. [S] see Eq. (|C9|I . The dot 
represents the xf^ vertex. 

In this formulation, the NLO contribution to the self- 
energy is given by the single one-loop diagram of Fig. [51 
corresponding to the expression: 



^NLO 



{x,x') = -G{x,x')D{x,x') 



(C2) 
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FIG. 9: The auxiliary field self-energy at LO in the 1/A'' ex- 
pansion. 



where G is the propagator of the field ip at LO and 

D{x,x')^{Tcx{x)x{x')) (C3) 

is the propagator of the composite field x- this order 
of approximation, the x self-energy is given by the dia- 
gram of Fig. [3J which corresponds to —^G'^{x^x') — 
^n(a;,x'), in terms of the one-loop bubble function 
n(a;,x'); see Eq. dTT]). One thus has 



iD-^{x,x') 



iN r 



6^°\x,x')-il].{x,x') 



(C4) 



where the first term on the right-hand side is the tree- 
level contribution; see Eq. (|C1|) . Using Eq. ([75]) . it is easy 
to check that the propagator D{x, x') is thus related to 
the infinite series of bubbles /(a;, x') as 



D{x,x') ^ ^[5'^''\x,x') + il{x,x') 



(C5) 



Inserting this expression in Eq. (jC2[) . one recovers the 
resuh ([751)- ([75]). 

Introducing the conformally rescaled correlator as 

D{x, x') = [a{'n)a{in')]-'^ D,{tj, t?', |X - X'|), (C6) 
one gets, in comoving momentum space, 

D,{r^,r^\K)^^[dc{v-r/)+tL{r^,f^',K)\ . (C7) 

The p-representation of the auxiliary field correlator is 
obtained as 

D,{rj,r^',K)^KD{p,p') (C8) 
with p = —Krj and p' = —Krj'. One has, in particular. 



D{p,p') 



iX_ 

3N 



'5Ap^p') + ii{p,p') 



(C9) 



Finally, we mention that the 2PI 1/A^ expansion, dis- 
cussed above, can also be formulated in the auxiliary field 
formalism .54, ^5] ■ The relevant equations are formally 
the same as presented here, with the understanding that 
the correlators G and D are the NLO ones — instead of 
the LO ones — to be determined self-consistently from the 
SD equations at NLO. 
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